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AERONAUTIC SYMBOLS 

1. FUNDAMENTAL AND DERIVED UNITS 


w 

Q 

m 

1 

v 

S 

S„ 

G 

b 

c 

A 

V 

2 
L 
D 
D> 

Dt 

D P 

G 



Symbol 

Metric ^ 

English 

Unit 

Abbrevia- 

tion 

Unit 

Abbrevia- 

tion 

Length. 

Time 

Force 

l 

t 

F 

meter 

second — _ 

weight of 1 kilogram _ _ 

m 

s 

kg 

foot (or mile) 

second (or hour) 

weight of 1 pound 

ft (or mi) 
sec (or hr) 
lb 

PnwPT 

P 

V 

horsepower (metric) 


horsepower. _ 

hp , 

mph 

fps 

1 U'Y VI — — — — — — — 

Speed 

(kilometers per hour 

\meters per second . 

kph 

mps 

miles per hour. 

feet per second _ — 


2. GENERAL SYMBOLS 


Weight=ra# 

Standard acceleration of gravity=9.80665 m/s 2 
or 32.1740 ft/sec 2 
W 

Mass=— 

g 

Moment of inertia= mP. (Indicate axis of 
radius of gyration k by proper subscript.) 
Coefficient of viscosity 


v Kinematic viscosity 

p Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-m""*-s 2 at 15° C 
and 760 mm; or 0.002378 lb-ft~ 4 sec 2 
Specific weight of “standard” air, 1.2255 kg/m 3 or 
0.07651 lb/cu ft 


3. AERODYNAMIC SYMBOLS 


Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio, -g- 
True air speed 
Dynamic pressure, ~ pV* 

Lift, absolute coefficient C L =-^ 


<L S 

D 


Drag, absolute coefficient C d—^ 
Profile drag, absolute coefficient 


__A 


D t 


Induced drag, absolute coefficient C Di = ^ 
Parasite drag, absolute coefficient C Dp = ^ 

Q 

Cross-wind force, absolute coefficient Cc— 


Angle of setting of wings (relative to thrust line) 
i t Angle of stabilizer setting (relative to thrust 
line) 

Q Resultant moment 

ft Resultant angular velocity 

R Reynolds number, p-^ where l is a linear dimen- 
sion (e.g., for an airfoil of 1.0 ft chord, 100 mph, 
standard pressure at 15° C, the corresponding 
Reynolds number is 935,400; or for an airfoil 
of 1.0 m chord, 100 mps, the corresponding 
Reynolds number is 6,865,000) 
a Angle of attack 

e Angle of downwash 

a 0 Angle of attack, infinite aspect ratio 

a< Angle of attack, induced 

a a Angle of attack, absolute (measured from zero- 

lift position) 

y Flight-path angle 
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BENDING AND SHEAR STRESSES DEVELOPED BY THE INSTANTANEOUS ARREST OF 

THE ROOT OF A MOVING CANTILEVER BEAM 

By Elbridge Z. Stowell, Edward B. Schwartz, and John C. Houbolt 


SUMMARY 

A theoretical and experimental investigation has been made 
oj the behavior of a cantilever beam in transverse motion when 
its root is suddenly brought to rest. Equations are given for 
determining the stresses, the deflections, and the accelerations 
that arise in the beam as a result of the impact. The theoretical 
equations, which have been confirmed experimentally, reveal 
that, at a given percentage of the distance from root to tip, the 
bending stresses for a particular mode are independent of the 
length of the beam whereas the shear stresses vary inversely 
with the length. 

INTRODUCTION 

When an airplane lands, the vertical component of the 
velocity is rapidly reduced to zero. In the absence of a 
thorough analysis of the stresses that arise from such shocks, 
it is customary for engineers to assume that the landing 
loads are static and independent of the elastic properties of 
the structure. As an initial step in the study of elastic 
structures under shock loads, an investigation has been made 
to determine the effect on a simple structure of the sudden 
arrest of its motion and the effect of the geometry of the 
structure on the stresses that result. The particular case 
treated in this report covers the basic problem of the in- 
stantaneous arrest of the root of a moving cantilever beam. 
The solution of this problem gives the energy consumed in 
exciting the different modes of vibration and the stresses, 
deflections, and accelerations that result throughout the 
beam. 

This investigation is based on the usual engineering 
beam theory in which the deflections are considered to be the 
result of bending alone and shear deflections are neglected. 
The theory, as applied to ordinary beams, gives reasonably 
good results as long as the distance between inflection 
points is greater than a few times the depth of the beam. 
When this theory for beam action is used in vibration 
problems, such as the problem in the present paper, the 
results are satisfactory for those modes of vibration for 
which the nodes are not too close together. This report 
summarizes the results of a theoretical solution, given in 
the appendix, and presents an experimental verification 
of these results. 

SYMBOLS 

E modulus of elasticity 

7 weight density of material 

X coefficient of equivalent viscous damping of 

material 

c velocity of sound in material 



g 

L 

I 

A 

p 

x 

V 
t 

V 
n 

K 


acceleration of gravity 
length of beam 

moment of inertia of cross section of beam about 
neutral axis 

cross-sectional area of beam 
radius of gyration of cross section of beam 



coordinate along beam measured from root 
distance from neutral axis of beam to any fiber 
time, zero at impact 


operator 



integers 1, 2, 3, etc., designating a particular 
mode of vibration 
nth positive root 1 +cos 6 cosh 0=0 
undamped natural angular frequency of nth mode, 

radians per second ( pc j n . 


damped natural angular frequency of nth mode, 

radians per second (^ n ^/ ^When 
\2 2 

the “frequency” is defined by 





X 2 “n 2 
4 E 1 


) 


v velocity of beam prior to impact 

w(x, t) deflection of beam at station x and time t 
w n (x, t) deflection of beam at station x and time t for 
nth mode of vibration 

a{x, t) acceleration of beam at station x and time t 
a n (x , t) acceleration of beam at station x and time t for 
nth mode of vibration 

a(x, y, t) bending stress in beam at station x , distance 
from neutral axis y, and time t 
y, t) bending stress in beam at station x, distance 
from neutral axis y, and time t for nth mode of 
vibration 

r(x , t) average shear stress over cross section of beam at 

station x and time t 

f n (x, t) average shear stress over cross section of beam 

at station x and time t for nth mode of 
vibration 

A n bending-stress coefficient 

B n shear-stress coefficient 

C n deflection coefficient 
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RESULTS ANI) CONCLUSIONS 

THEORETICAL 

When a cantilever beam under uniform translation in 
a direction perpendicular to its length lias it root instanta- 
neously brought to rest ; there is excited a theoretically infi- 
nite number of modes of vibration. With each successive 
mode, damping has an increasing influence upon the fre- 
quencies and amplitudes of vibration and, for sufficiently high 
modes, even changes the type of motion from oscillatory to 
nonoscillatory motion. In the lower modes, however, damping 
has little effect, and only terms of the first order in damp- 
ing need to be included in the equations. Only the equa- 
tions applicable to the lower modes, which alone are of 
importance in any practical ease, are presented in this section 
of the paper. For a more complete treatment of damping, 
see the appendix. 

The angular frequencies (27r times the frequencies in cps) 
are given by the equation 

( 1 ) 

where 6 n has the following values for successive modes of 
vibration: 

01=1.875104 05 = 14.137168 

0 2 =4.694O98 06 = 17.278759 

03=7.854757 \ f 

04=10.995541 6n ~2 ( 2n_ l )*> n > Q 

The energy that the beam possesses before impact is 
consumed in exciting the various modes of vibration and 
is distributed among the modes as follows: 


Mode, n 

Percentage of 
energy 

1 

61.3 

2 

18.8 

3 

7.4 

4 

3.2 

5 

1.9 

6 

1.3 

7 to 00 

6.1 


This distribution of energy among the different modes of 
vibration is presented graphically in figure 1. 



All stresses, deflections, and accelerations are damped 
sinusoidal functions of time and vary along the length of 
the beam. The bending stress a n (x, y , t) and the average 
shear stress r n (x, f), associated with the nth mode of vibra- 
tion, are given by the equations 

Xoj„- . 

V 2E 1 

2/j t)^=A n ~~ Ec sin (2) 

c p 

Aotfn 2 , 

2E 1 

t „(x, t) = B n ~j^ Ee sin a >„t (3) 


The variation of the dimensionless coefficients A„ and B„ 
with ar/L is given for n— 1, 2, and 3 in figures 2 and 3. The 



highest values of A n and B n , and hence the highest stresses, 
occur at the root of the beam. These values, for the first 
six modes, are 


. Mode, n 

A n at root 

J3 n at root 

1 

1. 566 

2. 146 

2 

.868 

4. 149 

3 

.509 

3. 994 

4 

. 364 

4. 00 

5 

. 283 

4.00 

6 

.231 

4. 00 


The foregoing values of A„ and B n at the root are presented 
graphically in figure 4. 

The maximum values with respect to time of a n (x, y, t) 
and r n (.r, t) associated with the nth mode of vibration, when 
the effects of damping are neglected, are 
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3 



a n {x,y)=A tt ^E (4) 

c p 

r n {x)=Bn- c J j E (5) 

The deflections w u (x y t) for the nth mode of vibration are 
given by the equation 

T2 Xco " 2 , 

. „ V U 2E • , i r»\ 

w n (x, t) — C n ~ — e sin c o n t (6) 

C p 

The accelerations a n (x, t) for the nth mode, when damping 
is sufficiently small, are given by 

a n (x, t)=~a) n 2 w n (x, t) (7) 

The variation of the dimensionless coefficient C n with xjL is 
given for n= 1, 2, and 3 in figure 5. 

The equations (4) to (7) for stress, deflection, and accelera- 
tion give the values associated with the nth mode of vibra- 
tion. Since all modes of vibration occur simultaneously, the 
net results are the superposition of the effects of all modes. 
This superposition gives the following equations: 



Mode 



Figure J. — Values of bending-stress coefficient An and shear-stress coefficient B n at root. 
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For bending stress, 

/ Xwi 2 , 

v y ( — 1 

a (z> y j 0 ~ E I A x e ' sin w x t 

+A 2 e 2E sin a) 2 t-\- . . . ) 

For average shear stress, 


/ Xa;|2 ^ 

t(z, t)= ~ j^EyBie 2E sin uj 


X oj2~ 

4“ 2E sin 4“ 


2E ^ ^ 


For deflection 


wOr, 




Xcji ? 

CV 2 # sin co^ 


Xo?2- _ \ 

-)-(7 2 e 2 £ sin tojf+ . . .) 


(S) 


(9) 


00 ) 


For acceleration, when damping is sufficiently small, 

a(x ’ t)= l 7 ( c “ 

Ams- \ 

sin co 2 £+ . . . J (11) 


_Xm,2 / 

\u\ 2 e 2 k sin 

Xojo- 


The equation for bending stress (equation (4)) reveals 
that, at a given percentage of the distance from root to tip, 
the bending stress for a particular mode is independent of 
the length of the beam and depends only on the velocity 
before impact. The equation for shear stress (equation (5)) 
reveals that the shear stresses at any station vary inversely 
with the length of the beam. These results are contrary to 
those that might be expected on the basis of experience with 
the static behavior of structures. For this reason an ex- 
perimental investigation was made. 


EXPERIMENTAL 

A circular steel tube of 1-inch outside diameter and 0. 028- 
inch wall thickness was mounted symmetrically on the end 
of a pendulum to form a pair of cantilever beams. (See 
fig. 6.) The pendulum was permitted to start its swing from 
a predetermined position and was suddenly brought to rest 
at the bottom of its swing against an electromagnet used to 
prevent rebound. The effect of length was studied by 
reducing the length of the tube in successive tests. The 
bending and shear strains were measured by electrical strain 
gages that were mounted on the tube as shown in figure 7. 
Each pair of gages was incorporated into a Wheatstone 
bridge circuit as shown diagrammatically in figure 8. The 
outputs of the bridge systems were fed through a strain-gage 



Figure t>.— Pendulum assembly used in impact test. 



amplifier into a multichannel oscillograph that recorded the 
strains on moving photographic paper. The amplitude of 
the components of strain due to the modes of higher fre- 
quency was reduced, however, because of the response charac- 
teristics of the oscillograph. The frequency-response curve 
for the oscillograph used is given in figure 9. 

Typical records for tubes of two lengths are shown in 
figure 10. Inspection of the record for the cantilever beam 
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O 400 600 1200 1600 2000 2400 


Frequency t cps 

Figure 9— Frequency response of strain-gage amplifier and Miller (2000~) oscillograph; 

0.69 critical damping. 

26% inches long shows the superposition of the second and 
third modes upon the first mode. The record shows that, in 
the case of the bending strain, the contribution of the second 
mode is small; whereas, in the case of the shear strain, the 
contribution of the second mode is large. This observation 
confirms qualitatively the theoretical results shown in figure 4. 
The same effect is not shown, however, in the record for 
the cantilever beam 11% inches long because of the combined 
action of damping and reduced response of the oscillograph 
to the higher frequencies associated with this short length of 
tube. 

The bending stresses computed by use of equation (8), in 
which only the first three modes are used, are given by the 
solid-line curve of figure 11 for the cantilever beam 26% inches 
long. Comparison of this curve with the record obtained 
during the first % cycle of the first mode (see fig. 10) shows 
good agreement as regards the wave shape. 

Because of the damping present in the tube and the 
response characteristics of the oscillograph, the only com- 
ponent of vibration that could be satisfactorily recorded for 
all lengths of cantilever tube was the fundamental or first 
mode. The quantitative results of the tests consequently 
were based upon this mode of vibration. This procedure is 
sound because the effects of the various harmonics are 
independent of one another. In the analysis of the results, 
the data had to be corrected for the influence of the magnet. 



(a) Cantilever length, 26H inches. 

(b) Cantilever length, 11 H inches. 


Figure 10. — Portions of typical records obtained for two different lengths of tube. 



Figure 11.— Theoretical wave form for extreme-fiber bending stress at root obtained from 
the first three modes of vibration. Steel tube, 1-inch outside diameter; wall thickness, 
0.028 inch; length, 26H inches. 
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(a) Bending. 

(b) Shear. 

Figure 12. — Comparison of experimentally and theoretically determined maximum stresses 
of fundamental mode at root of cantilever tube. Impact velocity, 2.5 feet per second. 


The observed frequencies are compared with the frequen- 
cies computed from equation (1) for the first mode in the 
following table: 


Length (in.) 

Frequency 

Observed 

(cps) 

Computed 

(cps) 

46$4 

17.5 

17.5 

::•••: 

27.9 

28.2 

2624 

52.1 

53.2 

1624 

131 

137 

1124 

272 

277 


The experimental values of extreme-fiber bending stresses 
and the shear stresses at the root, for the fundamental mode, 
are plotted in figure 12. In figure 12 are also shown the 
corresponding theoretical curves of equation (4) for bending 
and equation (5) for shear with n taken as 1. It is ob- 
served that the experimental points follow the trend of and 
lie close to the theoretical curves. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., September 27, 1944 • 


APPENDIX 


THEORETICAL DERIVATION 


General analysis. — Consider a beam of uniform cross sec- 
tion in equilibrium. If a portion of the beam is suddenly 
disturbed, as by a shock, in a direction perpendicular to its 
length, the beam is set into damped bending oscillations. 
The equation of motion for these bending oscillations is given 
by the differential equation (reference 1) 


, „ d 4 w d 5 w 7 d 2 w 

hp d.ri +Xp 'dr'dl + ! 7 dt> =0 


(Al) 


rhe damping term Ap 2 is derived on the assumption 


that the longitudinal damping force per unit area at any 
point on the cross section of the beam is proportional to the 
rate of change of longitudinal strain at that point. (See 
reference 2.) This type of force is analogous to ordinary 
viscous drag, in which the tangential force per unit area is 
proportional to the rate of change of shear strain. With the 
Eo 

use of the notation c 2 =^> equation (Al) can be written 


d 4 w , X d 5 w . 1_ d 2 w 

dt 2 


(A2) 


In accordance with the Heaviside operational methods 
(reference 3), equation (A2) may be reduced to an ordinary 

differential equation of the fourth order by writing 


thus, 


/ 1 , A \d*w_.p* 


c 2 p 2 


w= 0 


(A3) 


The general solution of equation (A3) is 


w=P cosh 6 J j- -\-Q sinh 6 j ■ +i? sin 6 cos 6 j (A4) 


where 

0=L I V — 
\ Wyl+P g 


The coefficients P, Q, R , and S are to be determined from 
the boundary conditions. The case under consideration is 
that of a cantilever moving with uniform velocity v and 
having its base brought instantaneously to rest. The 
boundary conditions for this case are 



The velocity of the root as given by the first boundary con- 
dition is represented graphically in figure 13(a). The rules 
of the Heaviside calculus, however, have been devised for 
a disturbance, called the unit function 2, shown in figure 
13(b). By the principle of superposition, the velocity func- 
tion shown in figure 13(a) may be considered as a super- 
position of those shown in figures 13(c) and 13(d). The 
velocity therefore consists of a constant velocity?? (fig. 13(c)) 
added to the solution of the problem obtained by the Heavi- 
side expansion theorem for the disturbance shown in figure 
13(d). On the basis of this procedure, the first boundarv 




v- vl 

0 

t 

vl 

0 



V 

0 




T 
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condition may be written 


(10 , =p (w) *- o=_ci 


' x=0 


With the application of the boundary conditions to equation (A4), the operational form of the solution for the velocity 
(that induced by the disturbance) is found to be 

pW= 2TlTcoThVco^)[ (1 + cos e cosh e) ( cosh d z +c 0 s 9 z) 

+ sin 6 sinh 6 ^cosh 6 j— cos 9 (sinh 6 cos 9 +cosh 0 sin 9) ^sin 9 sinh 9 
Interpretation of this operational expression and addition of the constant velocity v gives for the total velocity 


bw(x } t) 
dt 


--v—vl-\-2v^2F (e n 


e 2 E 1 


COS 03 n 't~ 


Xc On 

2 E 


where 

e n 

- r 6 ” 2 

03 n — pC jr 2 

W,'=«, -y/ 


/ , X V 

V 


sin co.'f 


nth positive root of 1 + cos 6 cosh 0 = 0 

undamped natural angular frequency of nth mode, radians/sec 


XV, 


1— 4 ^-,” damped natural angular frequency of nth mode, radians/sec 

( x x\ / x x\ 

cosh 9 n j— cos 9 n j— (cosh 9 n sin 9 n + sinh 9 n cos 9 n ) (sinli 9 n j — sin0 n -^ j 

\ Un L)~~ 9 n (cosh 9 n sin 9 n — sinh 9 n cos 9 n ) 


(AS) 


(A6) 


Integration of equation (A6) with respect to the time with 
the condition (w ) t=0 = 0 gives for the deflection 


w(x, t) = 2v X 

n = 1 


X 6n i) -w' • ,, 

— — sm co n 7 1 


=- c u p ^ n 


V- 


where 


X 2 co tt ! 
' 4£ 2 


’(*•!) 


e 2£: sinw/^2 (A7) 


2 F 

Cn = ~a~2 

v n 


The contribution of the nth mode to the deflection is 

XcOn^ 

W„(x, 0=7 L c n — == == = e lh sin w „'tl (A8) 
c P I, \~u» 


* /,_xv 

V 1 4 E~ 2 


Xw 


When equation (A8) may be put in the form 


W n (Xy f) 

C P 


/ xw 

V 4£ 2 ' 


e 2, ' ; sinh u n 't 1 (A9) 


where now 

, /X 2 co „ 2 

"« — 1 w «y 4 jrv’2 ~ 1 

The form indicated by equation (A8), where <1, is 

characteristic of the lower modes and represents damped 
oscillatory motion. The form indicated by equation ( A9), 

where > 1 (damping greater than critical), is characteristic 

of the higher modes and represents subsidence motion. 

From equation (A6) for velocity and equation (A7) for 
deflection, the complete behavior of the cantilever may be 
determined. The quantities of interest are the bending 
stresses, the shear stresses, and to some extent the accelera- 
tions. When damping is present, the equations representing 
the contribution of the nth mode to these quantities may be 
given in the two forms indicated by equations (A8) and (A9). 
In subsequent equations, however, only the form indicated 
by equation (A8) is given because it is characteristic of the 
modes that are of practical importance. 

Bending stresses. — The bending stresses <j(x,y,t ) at any 
fiber distance y from the neutral axis are 

, ., ^ b ! w 

a (x, y, t) = Ey ^ 


-rr*v + a 

- E ~c P £ - 4 - 


1 


/ , xV 2 
V 4 E~ 2 


e “ " sin w n 't 1 


where 


A n = 2 
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/ • • Z' • x • x^\ 

sin 0 n sinli 6 n (cosh 6 n jr+cos 0 V j )— (cosh 9 n sin 0 n +sinh0 w cos 0 n ) (sinh 6 t] j+sin 0 n ^ ) 

6 n (cosh d n sin 0 n ~ sinli 0 n cos 0 n ) 


The bending stress due to only the nth mode is 


0"n (*G V) 0 — B r A n 

L p 


1 


Xo»n 2 . 
~~2E C 


V 1 - 


X'co* 2 

4E 2 


sin u n 't 1 


Shear stresses. — The average shear stress over the cross section r (x, t) is 


f ( x , t ) =Ep 2 


d 3 w 

dF 


0 -Li 72 = 1 


V 1 


xw 

' 4E 2 


Xcj,, 2 

„ 2E • r 

0 sm o) w t 1 


where 


sin 


B n = 2 


sinli ^sinh 


— sin 0 n jV (cosh 0 n sin 0 n +sinh 0* cos 0 W ) ^cosli 0 W ^+cos 0 n 
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The average shear stress due to only the nth mode is 
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Accelerations. — From equation (A6), with the aid of the relation 

pF{t) l = F{0)pl+F'{t) 1 


the acceleration anywhere on the beam is found to be 
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With the aid of the orthogonal properties of the functions F 6 n it is possible to show that the quantity 2 (®n 1 

reduces to zero when 0<|gl. At f=0, the quantity if^F (d„ j) equals zero, and only the term —vpl remains. 
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This term indicates that at f=0 an infinite acceleration of zero duration exists at the root. 

The acceleration due to only the nth mode is 
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Comparison with the expression for w n (x,t) (equation (A8)) 
shows that the acceleration for each mode is out of phase 
with the deflection. When damping is sufficiently small, 
however, the relation between the acceleration and the 
deflection reduces to the well-known result for undamped 
vibration 

Q / n(%) 0 ~ Wj^Wjt^X, t) 
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Positive directions of axes and angles (forces and moments) are shown by arrows 
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4. PROPELLER SYMBOLS 
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Power, absolute coefficient Cp =— Jfp 

5 1 pV* 

Speed-power coefficient =Jp^ 2 
Efficiency 

Revolutions per second, rps 
Effective helix angle = tan" 


5. NUMERICAL RELATIONS 


1 hp= 76.04 kg-m/s=550 ft-lb/sec 
1 metric horsepower=0.9863 hp 
1 mph=0.4470 mps 
1 mps=2.2369 mph 


1 lb=0.4536 kg 
1 kg=2.2046 lb 
1 mi= 1,609.35 m=5,280 ft 
1 m=3.2808 ft 



